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ABSTRACT 

Asteroid 21 Lutetia, visited by the Rosetta spacecraft, plays a crucial role in the reconstruction 
of primordial phases of planetary objects. Its high bulk density and its primitive chondritic crust 
( Weiss et al.||2011 ) suggest that Lutetia could be partially differentiated. 

We developed a numerical code, also used for studying the geophysical history of Vesta 
( Formisano et al., [submitted ), to explore several scenarios of internal evolution of Lutetia, differ- 
ing in the strength of radiogenic sources and in the global post-sintering porosity. 

The only significant heat source for partial differentiation is represented by ■^^Al, the other 
possible sources (^"Fe, accretion and differentiation) being negligible. 

In scenarios in which Lutetia completed its accretion in less than 0.7 Ma from injection of 
^^Al in Solar Nebula and for post-sintering values of macroporosity not exceeding 30 vol.%, the 
asteroid experienced only partial differentiation. 

The formation of the proto-core, a structure enriched in metals and also containing pristine 
silicates, requires from 1 to 4 Ma: the size of the proto-core varies from 6 to 30 km. 



Subject headings: planetary systems - minor planets, asteroids: 
satellites: formation - planets and satellites: interiors 



individual(21 Lutetia) - planets and 



1. Introduction 

Asteroid 21 Lutetia plays an important role in 
the comprehension of the origin and evolution of 
planetary objects. As pointed out by [Bottke et| 



nario in which 21 Lutetia experienced a partial dif- 
ferentiation with the formation of a metallic "core" 



al. (2005), the size of Lutetia allows its survival 



against impact disruption and that is why it pre- 
serves its original large scale structure. Data pro- 
vided by the Rosetta spacecraft s uggest an high 
bulk density (3400±300 kg m'^) ( [Patzold et af 
2011): this evidence, combined to the primitive 



overlain by a primitive chondritic crust (Weiss et 
ar][20lT| ). 

Geophysical and thermophysical history of 
Lutetia depends strongly on its initial composi- 
tion and global macroporosity, which are greatly 
uncertain. Considerations on the global structure 
lead some authors to consider Lutetia as those 



nature of the crust (carbonaceous or enstatite 
chondrites) (Coradini et al. 2011), depicts a sce- 



asteroids with abundant fractures and joints ( As- 
phaug|[2"009 ) for which the inferred macroporosity 
is in the range ~ 6 — 40% (Consolmagno 2008[ 



Wilkinson et al. 2002). Furthermore, all other 
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asteroids having similar size (except 20 Massalia) 
are thought to have macroporosities of > 5 — 10% 
and ranging up to ~ 80% ( Consolmagno|20()8 ) . A 
strict upper limit on Lutetia's macroporosity (i.e., 
^ 52%) is provided by a model assuming that 



the entire asteroid is below a very thin chondritic 
surface layer and made up of pure iron. 



Weiss et al. (2011 1 proposed three different sce- 



narios for producing the high bulk density of Lute- 
tia via partial differentiation. In the first sce- 
nario, primordial Lutetia has the same size as 
present and is initially undifferentiated. After ra- 
diogenic heating and subsequent internal melting, 
the metallic core and silicate mantle form with a 
decrease of the macroporosity and an increase of 
the bulk density. In the second scenario, primor- 
dial Lutetia has a larger radius and only a smaller 
volume fraction experiences the melting. The un- 
differentiated outer layer is removed by subsequent 
impacts while the bulk density increases. In the 
third scenario, the differentiation of Lutetia oc- 
curs early and at first the chondritic crust is not 
retained or it is removed by subsequent impacts: 
after this phase, chondritic debris are deposited on 
the asteroid. 

In order to investigate which of these scenarios 
is the most plausible, we studied the thermal his- 
tory of Lutetia. Before describing our model, how- 
ever, we want to review the works already present 
in literature, as a wide range of thermal models 
of planetesimals assuming ^^Al as the main heat 
source exists has been developed in recent years. 



Merk et al. ( 2002 ) analyzed the dependence of 



accretion rate on thermal evolution of planetesi- 
mal, concluding that accretion process had to be 
considered as long as the accretion time was not 
negligible respect to the half-life of the radionu- 
clides providing energy. The authors also made 
use of Stefan law formulation in order to incor- 
porate the effect of latent heat into thermal evo- 
lution, but they neglected the role of sintering 
and convection and they did not used a radiation 
boundary condition (i.e., they assumed a fixed sur- 
face temperature). 

Ghosh et al. (|2003 ) also focused their atten- 



tion on the importance of the accretion process 
in the thermal evolution of asteroids, by study- 
ing the case of 6 Hebe, showing the differences 
between instantaneous and incremental accretion 
cases. They used the model developed byjGhosF^ 
McSween (1998) with a moving boundary condi- 



tion and a radiation boundary condition. The au- 
thors concluded that incremental models provide 
a link between theoretical models of measurable 
quantities (i.e., peak temperatures, cooling rate. 



radioisotope closure times) in meteorites that were 
determined by their thermal histories. 



Hevey & Sanders (2006) incorporated convec 



tion in their thermal evolution models when the 
degree of partial melting exceeded 50 vol.%. They 
analyzed the effects of sintering, whose onset was 
set at about 700 K, starting with a high poros- 
ity (50%) and a low thermal conductivity. They 
considered instantaneous accretion and, as well as 
Merk et al. ( 2002 1 , they used a fixed boundary 



temperature. The main result of their work was to 
constrain the accretion of the parent bodies of dif- 
ferentiated meteorites to within the first 1.5 Ma, 
or at most 2 Ma, from CAIs. These planetesi- 
mals therefore formed before most chondritic par- 
ent bodies, which accreted at a later time. 



Moskovitz & Gaidos (2011 ) studied the thermal 



consequences of melt migration and, in particu- 
lar, they investigate how the redistribution of ^^ Al 
from the interior into a crustal layer would affect 
the thermal evolution of planetesimals. They con- 
sidered the case of instantaneous accretion exclud- 
ing sintering, convection and a radiation boundary 
condition, and they concluded that differentiation 
would be most likely for planetesimals larger than 
20 km in diameter and that accreted within ap- 
proximately 2.7 Ma from CAIs. 



Sahijpal & Gupta (2011) performed numerical 



simulations of the processes involving both aque- 
ous alteration and planetary differentiation. They 
used asteroids of 100 and 270 km as representative 
cases to study, starting with high porosity and low 
thermal conductivity. They inserted the sintering 
and simulated the convection (at 50% of melting 
of silicates) by raising the thermal conductivity by 
three orders of magnitude. Their model assumes 
a fixed boundary temperature and does not take 
into account Darcy's law formulation for the segre- 
gation of silicates and melt. They concluded that 
a convective molten iron core is necessary to ex- 
plain the remnant magnetism of the carbonaceous 
chondrite. 



Finally, Elkins-Tanton et al. (2011 1 investigated 



the possibility that early radiogenic heating of 
planetesimals could create partially differentiated 
bodies with a primitive crust and magnetic core 
dynamos. The undifferentiated crust must be 
thick enough to prevent the majority of impacts 
from breaching it and persist through the internal 
magma ocean stage. They concluded that plan- 



etesimals accreting before about 1.5 Ma after CAIs 
are likely to fully differentiate through radiogenic 
heating, while bodies that accrete past about 1.5 
Ma from CAIs would probably be characterized by 
a thick undifferentiated crust overlying a differen- 
tiated interior. 

In the present work we apply the thermal code 
we developed to study the thermal history of Vesta 
( jFormisano et al., submitted) in order to con- 
strain the formation time, the size and the mass 
of the proto-core of Lutetia by varying the time- 
delay in the injection of ^®A1 in the Solar Nebula 



(a similar approach is also presented in Hevey & 



Sanders (2006), Moskovitz & Gaidos (2011) and 



Elkins-Tanton et al. (2011 )) and the post-sintering 



macroporosity. We define as proto-core a structure 
enriched in metals, formed as a consequence of the 
metal percolation, but still containing pristine sil- 
icates. 



As in Hevey & Sanders (20061 and Elkins- 



Tanton et al. (2011), we assume instantaneous 



accretion and solve the heat equation with ra- 
diogenic heating provided by ^^Al and ^°Fe. As 
showed by Weidenschilling (2008), the formation 



of a planetesimal having the same size as Lute- 
tia requires about 10^ year or less. Moreover, as 



pointed out by Moskovitz & Gaidos (2011 ), recent 



dynamical studies treating the turbulent concen- 
tration of small particles in proto-planetary disks 



(Johansen et al. 2007 2009 Cuzzi et al. 2010) 



show that planetesimal can grow "nearly instanta- 
neously" in less than 100 yr to sizes of 100 km or 
larger. As a consequence, instantaneous accretion 
can be a reasonable approximation. 



Similarly to Moskovitz & Gaidos (2011), we 



do not analyze the sintering phase due to the 
large uncertainties associated to the assumed ini- 
tial porosities and to the corresponding thermal 
conductivities. 

We include the treatment of the latent heat 
based on the the Stefan law formulation, as in 



Merk et al. ( 2002 ), in order to incorporate its effect 



on the evolution and we use Darcy's law model for 
metal percolation. 

The developed scenarios will be useful to de- 
pict a reliable geophysical and thermal history 
and to constrain the internal structure of Lutetia. 
They also offer a complementary approach to the 
works present in literature, like the evidence for 



differentiation provided by gravitational features 
( Vincent et al.|20fl ) or by hydrocode modeling of 
the largest impact crater on Lutetia (lOklay et al.' 



2012) 



2. The Model 

2.1. Initial and Boundary Condition 

We consider 21 Lutetia as a spherical body of 
radius fixed to 50 km with initial homogeneous 
composition. We suppose that the object contains 
a metallic component (Y) (about 25 vol.%) and 
a silicatic one {X) (about 75 vol.%). This com- 
position is similar to those of H and L class of 
ordinary chondrites. The post-sintering porosity 
ranges from 10 to 30 vol.%. The initial tempera- 
ture (Tq), which is also the temperature of the So- 
lar Nebula in the surrounding of Lutetia, is fixed 
to 200 K ( |Lewis||1974"| . In Tab [I] we report the 
physical parameters used in our model. We use a 
radiation boundary condition at the surface and a 
Neumann boundary condition at the center (heat 
flux equal to zero): 



Tir,t = 0) = To 

r^i -0 

ldr\r=0 ~ ^ 
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where Tsurf is the temperature of the surface, e is 
the emissivity, K is the thermal conductivity, and 
cr is the Stefan-Boltzmann constant (see TabfTl). 

2.2. Physical and Numerical Description 

We numerically study the thermal evolution of 
21 Lutetia, solving the following coupled equa- 
tions: 



(pc), 



dT 
'~dt 



= V • A',„VT 



H, 



and 



dY_ 



VV = 0. 



(2) 



(3) 



Eq.ipl) is the heat equation with the radiogenic 
heat source (H) (Castillo- Rogez et al. 2007 ) , where 
(pc)m and Km are the overall heat capacity and 
the overall thermal conductivity, respectively. 
Eq.(p| is the advection equation which controls 
the metal percolation in the silicatic matrix, with 



the migration velocity being a function of the per- 
meabihty of the sihcatic medium, of the density 
contrast between molten metals and solid silicates, 
of the gravitational acceleration and of the viscos- 
ity of the molten iron ( Yoshino et al.|2d04 ). As in 
Ghosh & McSweenl ( 1998 ) , when melting tempera- 



ture of Fe-FeS (1213 K) is reached, the percolation 
of the metals takes place through the silicatic ma- 
trix and the formation of the proto-core occurs. 
To take into account in a simple way the latent 
heat during phase transition, the specific heat is 
modified through the Stefan coefficient: 



Ste 
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where x is the degree of melting: 



T-T. 



X = 



sol 



T, 



liq 
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(4) 



(5) 



and so: 

c^c{l + Ste), (6) 

assuming a linear growth of x with increasing tem- 
perature. 

The numerical method is based on ID, for- 
ward time-central space (FTCS), finite difference 
scheme in radial direction with Lax correction. A 
spatial grid (Ar) of 500 m is used. To ensure 
the stability of our numerical approach, we use an 
adaptive time increment according to the Courant- 
Friedrichs-Lewy stability conditions for each of the 
physical processes (heat diffusion, metal percola- 
tion, radiation boundary condition) we consider 



in our work. Following Toksoz & Solomon (19731, 



thermal conduction imposes the following critical 
time step: 



At, 



cond 






(7) 



In analogy with 'Toksoz & Solomon ( 1973 1, we can 



define the following critical time step associated to 
the radiation boundary condition: 



Atr 






(8) 



Finally, while the percolation of metals is taking 
place, we need to solve also eq.® and introduce 
a third critical time step: 



lALpQJ^Q 



Ar 



(9) 



where v is the velocity of the metal percolation. 
The Courant-Friedrichs-Lewy stability condition 
requires that the time step used in our model sat- 
isfies the following criterion: 



At < rain {Atcond, Atrad, Atperc) • 



(10) 



Therefore, at each temporal iteration of the pro- 
gram we select the actual time step basing on the 
minimum critical time step among those we com- 
puted. As the stability condition requires the ac- 
tual time step to be lower than the critical one, 
we chose to use a value equal to 90% of the small- 
est critical time step as a compromise between the 
competing needs for stability and performances. 
So, our time step is defined as: 

At = 0.9 X min {Atcond-, Atrad, Atperc) ■ (H) 

3. Results 

We explored several scenarios characterized by 
different strength of the energy sources (the ra- 
diogenic heat due to the decay of ^®A1) and val- 
ues of post-sintering macroporosity (10, 20 and 30 
vol.%). The scenarios are labeled NO (instanta- 
neous accretion, Atd = 0), Nl {Atd — 0.3 Ma) and 
N2 {Atd — 0.7 Ma). The main results are showed 
in TabJ2J in which we report the size, the time of 
formation and the mass of the proto-core (i.e., a 
structure enriched in metals and containing pris- 
tine silicates) , and maximum temperature reached 
after 5 Ma. In all the cases we analyzed, the 
maximum degree of sificate melting (~10 vol.%. 



corresponding to about 1450 K, Taylor] (1992)) is 



reached only in a limited region of Lutetia. Since 
the degree of silicate melting is very low and, fol- 



lowing Hevey & Sanders (2006), the onset of con 



vection should require that silicate melt fraction 
exceeds 50%, we conclude that heat transport via 
convection is negligible. In fact, even in those re- 
gions where a significant temperature difference is 
present, the Rayleigh number does not overcome 
the critical values lying between 1000 and 2000. 

In the case of post-sintering porosity of 10 
vol.%, we observe that in all scenarios Lutetia does 
not completely differentiate and only a proto-core 
forms (see the maximum temperature versus time 
profile of Fig[I];f)). In Fig[lja), after 0.1 Ma, in 
the three scenarios the temperature is lower than 
the solidus temperature of silicates and the aster- 
oid is homogeneously heated. After 0.5 Ma (see 



FigjlTb)), NO enters in the melting temperature of 
silicates, while in Nl and N2 the temperature are 
lower than 1425 K. After 1 Ma (see FigQc)), in 
Nl the temperature reaches the solidus tempera- 
ture of silicates, while in N2 the values are still 
low. In Figjljd), after 3 Ma, we observe a slight 
general increase of the temperature in all scenar- 
ios and then, after 5 Ma (see FigfTle)), in NO the 
temperature has almost the same value while in 
Nl and N2 we observe a general decrease. The 
proto-core size ranges from 25 to 31 km: the time 
of formation ranges from 1.2 to 2.2 Ma. 

The general trend, for a porosity of 20 vol.%, is 
similar to the previous case, but the values of tem- 
perature reached are lower because of the lower 
amount of material, and therefore energy sources, 
per unity volume (see FigJ2ra)). As we can observe 
in Figl2{b), after 0.5 Ma, the temperatures over- 
come the liquidus temperature of metallic compo- 
nent in NO and Nl, while in N2 the temperature 
is in the window of melting of metals. After 1.5 
Ma (see Fig|2jc)), Lutetia is in the heating phase 
for all the scenarios and, after 3 Ma (see FigMd)), 
the general trend is the same. In Figl2je), after 5 
Ma, we observe a general decrease of the temper- 
atures for Nl and N2. In Figl2jf) the maximum 
temperature versus time profile is reported. We 
observe that the maximum temperature is reached 
in the hottest scenario, characterized by instanta- 
neous accretion (NO), and the time of formation 
of the proto-core ranges from 1.7 to 3.2 Ma while 
the proto-core size ranges from 14 to 21 km. 

If we choose a value of porosity of 30 vol.%, 
we can observe that, after a isothermal phase (see 
FigMa)) for all scenarios, at 0.5 Ma the temper- 
ature overcomes the liquidus melting temperature 
of metals in NO, while in Nl it is in the windows 
of melting of metals and, in N2, it is lower than 
1213 K (see Figjgb)). After 1.5 Ma (see Figjgc)), 
Nl overcomes the liquidus temperature of metals 
and N2 enters in the windows of melting of met- 
als. The general trend is the same after 3 and 5 Ma 
(see FigMd) and (e), respectively). In no scenario 
the temperature reaches the solidus temperature 
of silicates (see FiglsFf)). The formation of the 
proto-core occurs from 2.3 to 4 Ma. Maps of Fig|4] 
summarize the results obtained for the three cases 
of porosity explored. 

In Fig|4[a-b-c) maximum temperatures are 
reached in the middle region of the asteroid (from 



30 to 40 km from the center) as a consequence 
of partial differentiation: in this region, in fact, 
there is less mass to be heated and so the tem- 
peratures are higher. In Fig|4Fd-e-f) we can see 
that the general trend is the same of FigQa-b-c) 
but the temperatures reach lower values than in 
the previous case because the porosity is increased 
(i.e., 20 vol.%). In NO the melting of silicates is 
possible, while in Nl and N2 only the melting of 
metals occurs. FigQg-h-i) show that high values 
of porosity (i.e., 30 vol.%) prevent the reaching of 
the melting temperature of silicates and in par- 
ticular, when the delay in the injection of ^^Al is 
greater (i.e., 0.72 Ma), the melting of metals is 
possible only in a narrow region of Lutetia, rang- 
ing from 5 to 15 km. 

Although we would expect an increase of tem- 
peratures with increasing porosity due to lowering 
of the thermal conductivity ( Opeil et al. 2010 1, the 



dominant effect of the increase in porosity is the 
decrease of the concentration of the radioactive 
source (^^Al) linked to the silicates: therefore, the 
hottest scenarios are those characterized by the 
lowest porosity, as we can see in Tabj2j 

We also observe that in all the scenarios ana- 
lyzed a surface layer ranging from 2 to 5 km is be- 
low the Curie temperature of carbonaceous chon- 
drites (corresponding to about 873 K): this means 
that a remanent magnetization of this body is pos- 
sible ( Richter et al.|[2012 |. The remanent magne- 
tization of Lutetia has not a certain nature: it 
could be externally generated by external sources 
in the primordial Solar System ( Weiss et al.|2010 ) 
or due to the presence of an internal core dynamo 
that could confirm the scenario of partial differen- 
tiation. 

We can make a direct comparison between our 
temperature profiles and those of other models. 
Comparing our results with |Hevey fc Sanders] 
( 2006 ) , in case of 50 km body and instantaneous or 
late accretion (0.75 Ma), maximum temperature 
reached (about 8000 and 4000 K, respectively) are 
very different from our values (in the hottest sce- 
nario of current paper, i.e., NO and 10 vol.% of 
porosity, we obtain 1455 K). Our thermal profiles 
also show a maximum in the "mantle" , due to the 
differentiation and the affinity of ^^Al with sili- 
cates. These global differences probably depend 
on the different methodology and initial conditions 
adopted in the two works. The authors concluded 



that the thin crust that formed is destined to be 
processed by the magma flowing beneath due to 
small impacts or convective drag. On the contrary, 



Elkins-Tanton et al. (2011) concluded that undif- 



ferentiated chondritic crust survive through the in- 
ternal magma ocean phase. If we define the forma- 
tion of magma ocean at 50 vol.% of silicate melt- 
ing, the small degrees of silicate melting reached in 
our scenarios prevent this situation and preserve 
the primitive unmelted crust. |Sahijpal fc Gupta] 
(20111 also showed the existence of a chondritic 



crust even but they assumed a fixed boundary 
temperature instead of a more realistic radiation 
boundary condition, so that the extent of the their 
chondritic crust was likely overestimated. Since 
the methodology assumed in |Moskovitz fc Gaidos 
(2011) is similar to that of our work, maximum 



values of temperature are more compatible with 



our owns than those of Hevey & Sanders ( 2006 1 , 



even if global profiles are different. 

4. Summary and Conclusions 

Observational data do not provide stringent 
constraints about the internal structure. Cur- 
rently we know that Lutetia possesses a chon- 
dritic crust (carbonaceous or enstatitic) and its 
high bulk density has been interpreted as in in- 
dication of the presence of a metallic core. The 
results of our model suggest that partial differen- 
tiation is possible: in fact, the maximum degree of 
silicate melting is about 10 vol.% in a limited re- 
gion of the "mantle" . This is consistent with the 
scenarios proposed by [Weiss et al. (2011), if the 
current macroporosity (10 - 30 vol.%) is the same 
as the post-sintering one. In all scenarios only the 
formation of a proto-core, a structure enriched in 
metals, occurs. The proto-core formation takes 
from 1 to about 4 Ma and proto-core size ranges 
from 6 to 30 km. The relative proto-core mass 
ranges from about 1 to about 36 % of the total 
mass. Our results suggest that the accretion time 
of Lutetia should not exceed 0.7 Ma from CAIs 
and the post-sintering macroporosity does not ex- 
ceed 30 vol.%. 

In all scenarios we considered, a primitive, un- 
differentiated crust survives thermal evolution 



(Elkins-Tanton et al. 2011 Sahijpal & Gupta 



2011 ) and it possibly could be reduced by sub- 



sequent impacts. 



A possible remanent magnetization is retained, 
if we assume a carbonaceous surface composition, 
since in all the scenarios we explored a surface 
layer of about 2-5 km is below the Curie temper- 
ature. 

The main source of energy is represented by 
^^Al while the contribution due to ^°Fe and other 
possible sources (e.g. accretion and differentiation 
processes) is negligible. 
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tance. M.F. thanks his friend Demetra De Cicco 
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Quantity 



Value 



Unit 



Reference 



Final primordial radius (i?) 

Density of metal (pmet) 

Density of silicate (psu) 

Specific heat of metal (solid) {cmet,soi) 

Specific heat of metal (liquid) {cmet.iiq) 

Specific heat of silicate(solid) (csu^soi) 

Specific heat of silicate (liquid) {cgnng) 

Latent heat of metal (L^net) 

Latent heat of silicate (Lsii) 

Metal solidus {T;;2f) 

Metal hquidus {TIJ"^*) 

Silicate solidus (T«^j) 

Silicate liquidus (Tf>') 

Thermal conductivity of metal (Kynet) 

Thermal conductivity of silicate (Kgn) 
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Table 1: Physical parameter values used in this work. 
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Table 2: Summary of scenarios. Size, time formation and mass of the proto-core and maximum temperature 
reached after 5 Ma are reported. 
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Fig. 1. — From top to bottom, from left to right, plots are: temperature vs distance from center at 0.1 Ma 
(a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs time 

(f), for (j)= 10 vol.%. 
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Fig. 2. — From top to bottom, from left to right, plots are: temperature vs distance from center at 0.1 Ma 
(a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs time 
(f), for 0= 20 vol.%. 
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Fig. 3. — From top to bottom, from left to right, plots are: temperature vs distance from center at 0.1 Ma 
(a), at 0.5 Ma (b), at 1.5 Ma (c), at 3.0 Ma (d), at 5.0 Ma (e) and maximum temperature profile vs time 
(f), for 0= 30 vol.%. 
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Fig. 4. — In the top line plots are: thermal history maps of NO (a), Nl (b) and N2 (c) for (/> = 10 vol.%; 
in the middle line plots are: thermal history maps of NO (d), Nl (e) and N2 (f) for (p — 20 vol.%; in the 
bottom line plots are: thermal history maps of NO (g), Nl (h) and N2 (i) for = 30 vol.%. Magenta isoline 
represents the Curie temperature of carbonaceous chondrites. 
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